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Abstract. 

It was predicted some time ago that the cavity dynamical Casimir effect (generation 
i-S^ , of photons from the initial vacuum state in a cavity with moving walls) might be 

observed if a boundary vibrates at the double frequency of some selected cavity mode. 
' However, to register the created photons one has to couple the cavity mode with some 

, detector. Considering the harmonic oscillator model of a detector, we analyze how 

^ , ' different coupling regimes can affect the statistics of the created quanta. 



PACS numbers: 42.50.Ar, 42.50.Lc, 42.50.Pq 



1. Introduction 



m 

^ ■ A possibility of creating quanta of the electromagnetic field from the initial vacuum 
' state in cavities with moving boundaries, nowadays called the Dynamical Casimir 
Effect (DCE), was a subject of numerous theoretical studies for a long time: see, 
e.g., the most recent reviews [TH3!- It was shown [H-E] that one might expect a 
^ ' considerable rate of photons generation inside ideal cavities with resonantly oscillating 
boundaries. The simplest model describing this effect takes into account a single 
resonant cavity mode whose frequency is rapidly modulated according to the harmonical 
law ojt = i^o[l + £sm{rit)] with a small modulation depth, \e\ ^ 1. We shall use 
dimensionless variables, setting h = uo = 1- Then the Hamiltonian for the resonance 
mode has the form [7] 

He = utn - ixt{a^ - a^'^), Xt = {^uJt)~^dujt/ dt, (1) 

where a and are the cavity annihilation and creation operators, and n = a^a is the 
photon number operator. It is well known that the number of photons created from 
the initial vacuum state is maximal if the modulation frequency is exactly twice the 
unperturbed mode frequency, i.e., rj = 2. The mean number of photons (n) and the 
Mandel factor Q = [{(Anf) - {n)]/{n) increase with time in this ideal (hereafter 
we use the subscript for the quantities related to the empty cavity) 

{no{t)) = smh\et/2), Q^{t) = 1 + 2{n^{t)). (2) 
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The field mode goes to the squeezed vacuum state with the following variances of the 
field quadrature operators x = {a + a'^)/\/2 and p = (a — a'^)/(\/2i) (in the system 
rotating with the frequency Uq = 1): 

1 1 
"^pp ~ 2^ ' '^^^ ~ 2^ ' ^^■^ 

But simple formulae ([2]) and (I3l) hold for the ideal empty cavity only. To register 
the emerging photons one has to couple the field mode to some detector. And here 
the problem of the back action of the detector on the field arises, because in many 
realistic cases the coupling between the field and detector can be much stronger than 
that between the field and vibrating cavity walls. This was noted in [8], where it was 
shown that for the simplest model of detector as a two-level 'atom', no photons can be 
created at all for the modulation frequency 77 = 2 if the field-atom coupling constant 
g is much bigger than the frequency modulation amplitude e. But the photons can be 
created if one adjusts the modulation frequency rj = 2(1 + n), choosing some nonzero 
(small) value of parameter k. 

Here we consider the model of the detector as a harmonic oscillator tuned to the 
same frequency as the selected field mode. Despite its simplicity, this model seems 
to be rather realistic in the case of the so called Motion Induced Radiation (MIR) 
experiment [9|[T0]. where the microwave quanta created via the DCE are supposed to be 
detected by means of a small antenna put inside the cavity. Since the inductive antenna 
(a wire loop) used in that experiment is a part of an LC-contour, it can be reasonably 
approximated as a harmonic oscillator. Therefore, the Hamiltonian describing the 
system under study (the field mode coupled to such an antenna) can be taken in the 
form 

H = a^a + b^b + g (ab^ + ba^^ - ixt [a^ - a^^) (4) 

where the coupling constant g is assumed to be real number. Of course, the quadratic 
Hamiltonian (jl]) is an approximation, since it does not take into account possible 
nonlinear phenomena, e.g., effects of saturation in the limit of very long times. Therefore 
it can be used under the condition eH <C 1. But in the present state-of-art experiments 
on DCE, the time scale ~ 1 (or slightly bigger) seems to be quite sufficient for our 
purposes. 

Hamiltonian (jlj) contains three real (small) parameters: g, e and k. Our goal is 
to find the domains in the space of these parameters where the photon generation is 
possible and to study different regimes of generation. Due to the interaction with the 
detector, the field mode appears in a mixed quantum state described by the statistical 
operator p. We are interested, in this paper, in the photon distribution function (PDF) 
/(m) = {m\p\m), where \m) means the mth Fock state of the field mode. For the initial 
vacuum states of the field mode and the detector, the time-dependent statistical operator 
is Gaussian. The general form of PDF of the Gaussian states is well known [TTHTG] . For 

I The term uJtn in ([T]) can be replaced by n because for |e| <C 1 the main effect of modulation is due 
to operators and a^^ in the squeezing part of He- 
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zero mean values of quadrature components x and p, it can be expressed in terms of the 
Legendre polynomials as follows [T7] : 

2^m/2 / 4A - 1 \ 

where 

L'± = l + 4A±2r, T = a^^ + app=l + 2{n), (6) 

/I \ ^ 2 

The Mandel parameter in the Gaussian states with zero first-order moments can be also 
expressed through the quantities A and (n) as 

Q = l + 2(n) -(A-l/4)/(n). (8) 

Another quantity we are interested in is the invariant squeezing coefficient p^] - [20] 



s = (9) 

which does not depend of possible rotations in the quadrature plane, being equal to 
unity for the vacuum or coherent states. 

2. Photon generation regimes 

The first two terms in Hamiltonian can be removed by going to the interaction 
picture. Besides, using the rotating wave approximation (RWA) we can remove rapidly 
oscillating terms in the product Xt {0? — a^^^ . Thus, we arrive at the new Hamiltonian 

h\^^^^ = -1(3 (a^e-^-* - at^e^-*) + g (afet + ha^) (10) 
with /3 = £:/4. The corresponding Heisenberg equations of motion 

da/dt = 2f3a^e^'''^ - igb, db/dt = -iga (11) 
can be solved analytically by means of the substitutions 

a(t) = e'^'d{t), b{t) = e'^'b{t), 
which result in equations with constant coefficients 

dd/ dt = 2l3a) — igb ~ ind, db/dt = —igd — iKb. (12) 



Looking for solutions to equations (1121) and their Hermitian conjugated partners in the 
form a, 6, a"!", fe"!" ~ e'^*, we arrive at the characteristic equation 

A4 + [k^ + g^- 2/32) + -g-'f- Ak^(3^ = (13) 
whose solution reads 
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The photon generation is impossible if Re(Ai 2,3,4) = for all four solutions ( IT^ . 
Otherwise, the real part of at least one characteristic value A is positive, meaning an 
exponential growth of solutions. Analyzing formula ( HM we conclude that the photon 
generation is impossible if the following three inequalities are satisfied simultaneously: 

K^ + g^> 2/32 (15) 

/3^-/3V + ^7V>0 (16) 

[K^-g'f >Ak'(3\ (17) 

If any of the inequalities (fT5l)-(lT71) is not satisfied, then an exponential growth of the 
mean number of photons can be observed. In figure [1] we show the regions in the 
parameter plane K-g where the photon generation from vacuum is possible. In this 
figure all parameters are normalized by /3 (i.e. formally we put /3 = 1). 




Figure 1. The regions (yellow in the online version) in the parameter plane K-g where 
the photon generation from vacuum is possible. 



3. Analysis of special cases 

3.1. Expected resonances for \k,\ = \g\ 

Condition (ITTj) is obviously broken if \k,\ = \g\, i.e. along the bisectrices in figured] 
The possibility of photon generation in this case seems quite natural, as soon as the 
corresponding modulation frequency rj = 2{1 + k) is exactly twice bigger than one of 
two eigenfrequencies u± = 1 ± g of the stationary part of Hamiltonian (jl]) (with = 0). 
Namely, this case was studied for the first time in [3j. In particular, for ^ \g\ 
photons can be created if \k, — g\ < \(3\. Under this condition the solutions to equations 
(fTTl) have rather simple explicit forms if, in addition, {Pt){P/g) ^ 1: 



a(^) = I [(ao - bo) cosh(/3t) + (aj - bl) sinh(/3t)] e'^' + ^ (ao + 60) e'''\ (18) 
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b{t) = - \{bo - ao) cosh(/3t) + (foj - aj) smh{/3t)] e'^' + - (oq + &o) e"*^* . (19) 



2 



2 



The mean numbers of quanta in both the modes coincide (for the initial vacuum states): 
Mt)) = {ubit)) = i smh\€t/A). (20) 

The photon generation rate in the field mode interacting with the oscillator detector 
turns out to be twice smaller than for the empty cavity [given by equation ([2])]. This 
result was obtained in [1] , but unfortunately the argument of the hyperbolic sine function 
there was twice bigger due to a misprint. The effect of diminishing the photon generation 
rate due to the resonance intermode interactions was discovered in [SlEI]- In the most 
strong form this effect manifests itself in effectively one- dimensional Fabry-Perot cavities 
with (quasi) equidistant spectra of eigenfrequencies [11122] • For other statistical properties 
of the field mode, we have the following formulae: 

A = ^ cosh2(/3t), Q = 1 cosh(2/3t) = 2{n) + ^, (21) 
(^xx = ^ cosh(/3t) [cosh(/3t) + sinh(/3t) cos(2^t)] , (22) 
(Tpp = ^ cosh(/3t) [cosh(/3t) - sinh(/3t) cos(2^t)] , (23) 



CTxp = 7^ cosh(/3t) sinh(/3t) sin(2^t). (24) 



1 
2 

The minimal value (with respect to fast oscillations with frequency 2g ^ f3) of the 
variance of any of two quadrature components is equal to 

Formula (E]) for the PDF in the field mode can be written in the case involved as (see 
also m) 

tanh(/3t) 



/(m) = itzrVT^^P^{-iz), z = , . (26) 

/4-tanh2(/3t) 



For big values of index m, one can use the asymptotical formula for the Legendre 
polynomials [23] 

P^(coshO ^ [-^^ ' /o([m + l/2]0, (27) 

where /o(2;) is the modified Bessel function. Taking into account known asymptotical 
formulae for the Bessel functions of big (complex) arguments and following the scheme 
described in |l7j, one can arrive at the formula 

[tanh(/3t)]'" [- tanh(/3t)]'" 

^(^^ ^ [2 - tanh(/3t)]-+i | [2 + tanh(/3t)]-+^ ^ ^^8) 



cosh(/9t) y vr {m + ^ 
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which is vahd under the condition m ^ 1 for both small and big values of the product 
Pt. It is worth comparing formula f l28|) with the strongly oscillating distribution 

in the squeezed vacuum state arising in the absence of interaction with the detector. The 
probabilities of observing odd numbers of quanta in the distribution (128) are close to zero 
if /3t <^ 1. But this case is not very interesting, since {n{t)) <^ 1 under this condition. 
In contrast, if /3t > 1, so that tanh(/9t) is close to unity and {n{t)) ^ exp(2/3t)/8, then 
one can rewrite ( l28l) as 

1 _|_ / 1 \mq— m— 1/2 

/(m) ^ ^ ^ - . (30) 

^27rm{n{t)) 

For m ^ 1 the second term in the numerator of fraction in formula ( 130|) is very small. 
Therefore this formula shows very smooth distribution, quite different from (!29|) . Note 
that for ^ 1 formula ( l29l) can be written (using the Stirling formula for the factorials) 
as fo{2k) ^ [7ik{n)]~^^'^. Comparing this expression with fl30l) for m = 2k we see that 
f{2k) ^ fo{2k)/2, so the distribution (l30l) can be considered as an average of even and 
odd values of the 'saw-tooth' distribution ( |29|) . The plots of exact distributions (!26|l 
and ( l29l) for m < 20 and (n) ~ 6, illustrating these observations, were given in [4]. 

On the basis of this example, one could suppose that the drastic change of the 
behavior of the PDF is due to the strong coupling with a detector, which plays a 
role of some 'reservoir' (note that thermal reservoirs usually cause 'smoothing' of any 
oscillatory behavior). However, the examples of the following subsections show that the 
real situation is more intricate, and even the strong coupling with a detector not always 
destroys the oscillations of the PDF or some other physical quantities. A rough analogy 
can be the case of nonthermal 'rigged' reservoirs, which can enhance oscillations of some 
functions. 



3.2. Surprising resonance at k = 

Figure [T] shows the existence of resonance photon generation for k = and for any value 
of the coupling constant g. This result, first discovered in [23], seems surprising, because 
in the absence of detector (for g = 0) the mean number of quanta in the case of a small 
detuning k 7^ is given by the following generalization of formula ([2]): 

{no{t)) = ^^J^2 si^h^ (V^V4 - K^) ■ (31) 

Formula f l3T|) shows that the deviation of the modulation frequency from the resonance 
value rj = 2 hj 2k, = e stops the photon generation in the empty cavity. Therefore 
it was natural to expect |3] that for g 0, the modulation frequency must be close 
to 2u}± = 2(1 ± g), with the deviation not exceeding something of the order of e. 
Nonetheless, in reality the photons can be created also for \k,\ < even if \g\ ^ \(3\. 
Perhaps, this happens due to some kind of quantum interference. The solutions of 
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equations of motion f llip with k = and arbitrary values of f3 and g can be found 
in 123] (similar equations were solved in the contexts of different other physical problems 
in |25l|26]). We bring here only some consequences of that solutions. The mean number 
of quanta in the field mode in the case of \g\ ^ \f3\ is equal to (for f3t 1) 



2/3^ 
1 + -sin(27t) + —-sm^ {-ft) 



(32) 



7 7 

where 7 = V — Again, the rate of photon generation is roughly twice smaller 
than in the empty cavity, but the mean photon number is approximately twice bigger 
than in the case of k, = g considered in the preceding subsection. Time dependences 
of the mean numbers of quanta in the field mode in different regimes are compared in 
figure El The third (blue) line from the left (corresponding to the case of k = 0) shows 
remarkable horizontal steps. This peculiar behavior was explained in 




Figure 2. The mean number of quanta in tlie resonance field mode versus the 
dimensionless time et for different parameters e, g and k. The order of lines from 
the left to the right: g = k = 0, e = 10^^; 13 = k — g — 10^^ (see subsection 13. 3|) : 
K ^ 0, g = lOe ^ 10-2; k = 5 = 10^2 = iQe. 
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Formula indicates that even-odd oscillations of the PDF can happen if the 
argument of the Legendre polynomial is close to zero [since -P2m+i(0) = while 
-P2m(0) = (— l)'^(2m — l)!!/(2m)!!]. In the most strong form these oscillations manifest 
themselves for pure quantum Gaussian (squeezed) states with 4A — 1 = 0, as one can 
see in formula (l29l) . In the case of k = we have 4A — 1 = 4(7^/3^ sin^(7t)/7^, and this 
quantity is very small if \g\ ^ \(3\. Therefore, contrary to the case of k = ^ /3, the 
PDF shows oscillations, and for m ^ 1 we have 



f(2r.) = tanh-"(/3^) (2m -1)!! 
' cosh(/3t) (2m)!! ' 

/(2m + 1) 2(2m + l)/32sin^(7t) 



/(2m) 



c/2sinh(2/3it) 



< 1. 



(33) 
(34) 



3.3. The intermediate regime k = g = P 

Another interesting special case admitting explicit solutions is (3 = k = g (the point 
of intersection of the bisectrice of the first quadrant and the circle + g"^ = 2/3^ 

[see condition ( [T5|) ] in figure [1]). The characteristic values are Ai,2 = ±v^/3 and 
A3,4 = ±i\/2/3, so that the following exact formulae hold (x = \/2 (5t): 

a{t) = —i=< oq Vs cosh(a;) — zsinh(a;) — isin(a;) 



X] 



+ al [3 sinh(x) + sin(x)] 



\/2 cos(a;) — \/2 cosh(x) — i sinh(a;) — i sin( 



6q sin(a;) — sinh(x) + iy/2 cosh(x) — i\/2 cos( 



X 



{a^a) = - 1 + 3 sinh^(a;) + sinh(x) sin(x) — cosh(a;) cos(x) 



(35) 



Despite the presence of trigonometric functions in formula ( 135|) . the mean number of 
photons grows practically exponentially without visible oscillations, as shown by the 
second line from the left in figure O The asymptotical rate of photon generation in this 
case is equal to e/\/2 - an intermediate value between e and e/2 characterizing the two 
adjacent curves. 

The parameters entering formula ([5]) for the PDF follows: 

4A — 1 = 2 cosh^(x) — 2 [sinh(x) sin(x) + cosh(x) cos(x)] 

1 2 
— - [cosh(x) sin(a;) — sinh(x) cos(x)] , 

= 8 cosh^(x) — 4cosh(x) cos(x) — - [cosh(x) sin(x) — sinh(x) cos(x)]^ , 



D 



For X ^ 1 we have 



-4sinh^(x) — 4sinh(x) sin(x) [cosh(x) sin(x) — sinh(x) cos(x)]" 



1 



4A- 1^-6^^(1-^/4), 



^ = (sin X — cos x)^. 
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^ 2e'%l - e/16), D_ ^ -e2-(l + e/8). 

Consequently, 



Since < ^ < 2, the argument of the Legendre polynomial varies from — z/yS to —i/Q, 
i.e. it cannot assume very small values. Therefore the PDF does not show noticeable 
oscillations, and can be well approximated (for 1 ^ m ~ (n)) by the formula pTflST] 

fM^'-^^i^^±mm. (36) 

The quantity showing oscillations in the case concerned is the invariant squeezing 
coefficient S = 2crmin- Indeed, since A ~ r ^ 1 for x ^ 1, formula can be 
simplified as 5 ~ 2A/r, so that S{x ^ 1) ~ (1 — ^/4)/3. Since ^(x) is a periodic 
function of time, the minimal quadrature variance amin does not go asymptotically to 
some limit, but it oscillates between the values 1/6 and 1/12. 



4. Conclusions 



The main results of this paper are as follows. We found the conditions of photon 
generation in a three-dimensional cavity with resonantly oscillating ideal walls when the 
resonance field mode is linearly coupled to a detector modeled as harmonic oscillator. 
The 'allowed' and 'forbidden' zones in the space of parameters e, g and k, are presented 
in figure [H We have shown that the main physical observables, such as the mean number 
of created quanta, their distribution function and the invariant squeezing coefficient, can 
show either smooth monotonic behavior or some kinds of oscillations, depending on the 
parameters characterizing the process. However, the oscillations of different quantities 
seem to be uncorrected, according to the examples considered. 
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